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FOREWORD 


Fast  Antiship  Cruise  Missiles  (ASCMs)  represent  one  of  the  most  severe  threats  to  surface 
ships  of  the  Navy.  The  radar  tracking  of  low-flying  ASCMs  is  made  difficult  by  the  multi- 
path  phenomenon.  The  random  nature  of  diffuse  reflections  from  the  sea  surface  and  internal 
thermal  noise  of  the  receiver  produce  errors  that  severely  degrade  the  captured  signal.  The 
diffuse  reflections  result  in  large  errors  in  the  monopulse  measurements  by  the  threat’s  ele¬ 
vation  angle.  This  report  describes  a  study  of  the  effects  of  the  various  sources  of  noise  on 
angle  measurements  by  monopulse  radar.  Quantifying  the  dependence  of  the  accuracy  of 
the  angle  measurements  on  the  noise  is  very  important.  The  monopulse  ratio  is  generally 
a  complex  quantity  because  of  arrivals  of  reflected  signals  from  different  directions.  Thus, 
the  bivariate  density  distribution  of  the  monopulse  measurements  along  both  the  real  and 
imaginary  axes  is  of  interest. 

The  present  study  was  supported  through  the  Office  of  Naval  research,  in  part  by  the 
Independent  Research  Program  of  the  Naval  Surface  Warfare  Center,  Dahlgren  Division, 
managed  by  H.  Raymond  Riedl  (R04),  and  in  part  by  the  Surface  Weapons  Systems  Tech¬ 
nology  Block  Program,  managed  at  NSWCDD  by  Robin  Staton.  This  technical  report  was 
reviewed  by  Dr.  Ronald  Helmick  (B32). 
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Systems  Research  and  Technology  Department 
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ABSTRACT 


The  accuracy  with  which  the  angular  position  of  an  object  can  be  measured  by  monopulse 
radar  depends  on  the  signal  level  relative  to  that  of  the  various  noise  sources.  The  signals 
captured  in  the  sum  and  difference  channels  of  a  monopulse  system  are  generally  represented 
by  complex  phasors  because  the  returning  signals  can  arrive  from  different  directions  with 
different  phases.  The  primary  variables  consisting  of  the  four  real  quantities  which  define 
these  phasors  are  stochastic  variables  because  of  the  noise.  Consequently,  the  Monopulse 
Ratio  (MR),  which  is  the  difference  phasor  divided  by  the  sum  phasor,  is  also  a  stochastic 
variable  because  it  depends  on  these  measured  quantities.  Gaussian  noise  in  the  sum  and 
difference  channels  does  not  usually  produce  a  MR  with  a  Gaussian  distribution. 

Straightforward  statistical  methods  are  used  to  obtain  the  joint  density  of  the  real  and 
imaginary  parts  of  the  MR  for  several  cases  under  the  assumption  that  the  primary  variables 
have  a  multivariate  Gaussian  distribution.  Three  cases  of  noise  correlation  in  the  sum  and 
difference  channels  are  discussed,  and  the  corresponding  joint  MR  densities  are  developed. 
These  represent  a  quantitative  description  of  the  effects  of  noise  occurring  in  the  primary 
channels.  The  marginal  density  representing  the  real  part  of  the  MR,  which  has  figured 
prominently  in  the  literature,  is  displayed  in  a  few  simple  cases.  The  results  of  Monte  Carlo 
experiments  are  given  to  show  consistency  with  analytically  determined  density  functions 
and  to  demonstrate  empirical  derivation  of  the  density  functions  in  cases  where  an  analytical 
approach  is  not  feasible. 

The  variance  of  the  MR  determined  by  a  single  observation  is  calculable  {i.e.,  finite)  only 
if  the  long  tails  of  the  density  function  are  somehow  limited.  Density  functions  with  finite 
variance  are  obtained  if  a  threshhold  is  placed  on  the  magnitude  of  the  sum  phasor  to  discard 
returns  lying  below  the  threshhold.  The  resulting  variance  is  developed  for  two  simple  cases 
where  the  primary  mean  values  are  zero.  The  dependence  of  the  joint  MR  density  on  a 
threshhold  value  is  studied. 
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CHAPTER  1 
INTRODUCTION 


An  amplitude-comparison  monopulse  radar  system  measures  target  elevation  (or  azimuth) 
by  using  two  antenna  lobes  which  are  squinted  (i.e.,  pointed  in  slightly  different  directions) 
relative  to  each  other.  The  system  emits  a  pulse  through  both  lobes,  and  combines  the 
returning  signal  received  in  the  two  lobes  into  sum  and  difference  channels.  A  target  situated 
on  boresight  (i.e.,  midway  between  the  two  lobe  axes)  will  reflect  identical  signals  into  the 
two  lobes  resulting  in  zero  output  in  the  difference  channel.  For  a  target  off  the  radar 
boresight,  the  output  in  the  difference  channel  will  not  be  zero.  The  output  of  the  sum 
channel  provides  a  reference  to  which  the  amplitude  and  phase  in  the  difference  channel  can 
be  related.  In  the  absence  of  glint,  multipath,  jammers  and  multiple  targets,  this  relative 
phase  will  be  zero  or  180  degrees  depending  on  which  side  of  boresight  the  target  lies.  The 
amplitude  of  the  difference-channel  voltage  relative  to  that  of  the  sum  channel  determines 
angular  distance  from  boresight  if  the  target  is  detected  within  the  radar  beam. 

In  practice,  the  relation  of  the  difference  signal  to  the  sum  signal  is  represented  by  means  of 
the  Monopulse  Ratio  (MR),  which  is  the  ratio  of  the  difference  phasor  to  the  sum  phasor.  In 
the  ideal  situation  described  in  the  previous  paragraph  the  MR  is  real,  but  in  actual  situations 
the  MR  is  complex  because  the  returning  signals  often  arrive  from  different  directions  with 
different  phases  within  the  main  beam  of  the  radar.  The  MRs  in  the  azimuth  and  elevation 
directions  are  converted  to  angles  in  these  directions,  often  considering  only  the  real  parts 
of  the  MRs. 

Several  situations  give  rise  to  the  directional  spread  of  the  arriving  signals.  The  most 
troubhng  and  severe  situation  arises  when  a  target  is  close  enough  to  the  sea  surface  that 
both  the  target  and  its  image  lie  within  the  main  radar  antenna  beam.  This  is  the  well- 
known  multipath  problem,  which  has  received  considerable  attention  over  the  years  without 
a  sufficiently  accurate  method  of  determining  target  elevation  having  been  found.  Multiple 
targets  and  jamming  also  produce  signals  arriving  from  different  directions.  Signals  also 
arrive  within  a  significantly-wide  cone  of  directions  from  a  single  extended  target  (e.g., 
glint),  even  in  the  absence  of  multipath.  All  these  phenomena  result  in  a  phase  difference 
between  the  sum  and  difference  channels  of  the  monopulse  radar  that  produces  a  MR  which 
is  complex. 

Originally,  these  phenomena  were  handled  by  determining  target  elevation  and  azimuth 
from  the  real  part  of  the  MR  only.  The  imaginary  part,  which  arises  because  of  the  non- 
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vanishing  phase  differences  of  the  arriving  signals,  was  considered  an  unavoidable  nuisance. 
However,  in  recent  years  the  imaginary  part  of  the  MR  has  been  used  to  obtain  more  accurate 
measurements  of  target  angle. 

The  primary  variables,  consisting  of  the  four  real  quantities  which  represent  the  phasors, 
are  stochastic  variables  because  of  random  errors  resulting  from  the  noise.  Consequently, 
the  MR,  which  depends  on  these  measured  quantities,  is  also  a  stochastic  variable.  It  is  of 
interest  to  consider  how  the  distribution  of  the  empirically-determined  MR  depends  on  the 
noise  sources.  The  noise,  particularly  that  component  resulting  from  diffuse  sea-reflection, 
is  difficult  to  quantify.  The  effects  of  noise  having  a  Gaussian  density  with  specified  means 
and  variances  are  studied.  This  assumption  permits  the  calculation  of  the  resulting  MR 
density  resulting  from  specific  variance  models.  In  view  of  the  uncertainties,  especially  in 
the  mechanism  of  reflection  from  the  rough  sea  surface,  it  does  not  seem  worthwhile  to  treat 
a  more  complete  model. 

The  real  and  imaginary  parts  of  the  monopulse  ratio  can  be  expressed  in  terms  of  the  real 
and  imaginary  parts  of  the  sum  and  difference  phasors  by  the  relation 


p-f  = 


X  +  iy 
u  -H  iv 


(1.1) 


where  p+iq  is  the  complex  MR,  x  +  iy  =  D  is  the  difference-channel  phasor,  and  u  +  iv  =  S  is 
the  sum-channel  phasor.  The  x  and  u  are  commonly  referred  to  as  the  in-phase  components, 
and  the  y  and  v  as  the  quadrature  components.  The  analysis  will  be  done  exclusively  in  the 
real  domain.  Accordingly,  the  relations 


xu  yv 


and 


yu  —  XV 


(1.2) 


define  the  stochastic  variables  p  and  q,  whose  distribution  is  sought,  while  the  distribution  of 
the  primary  quantities  x,  y,  u,  v  will  be  postulated  to  have  a  multivariate  normal  distribution 
with  specified  means  and  covariance  matrix.  Accordingly,  the  x,y,u,v  variables  are  allowed 
to  have  nonzero  means  and  a  general  covariance  matrix. 

Since  the  mapping  is  four  (be.,  primary  real  quantities)  onto  two  (be.,  secondary  real 
quantities),  a  multiplicity  of  density  functions  of  the  (x,y,u,u)  variables  will  result  in  the 
same  density  function  for  (p,  9).  The  ratio  (1.1)  is  not  affected  if  both  numerator  and 
denominator  are  scaled  by  the  same  complex  quantity.  Accordingly,  the  MR  is  determined 
solely  by  the  magnitude  ratio  and  phase  difference  of  D  relative  to  S.  The  formulation  is 
sometimes  represented  as  a  dependence  of  the  MR  on  three  primary  quantities,  which  might 
be  chosen  to  be  the  D  to  S  relative  phase  and  the  two  magnitudes.  It  is  usually  not  convenient 
to  reduce  the  number  of  primary  variables  further  by  replacing  the  two  magnitudes  by  their 
ratio  because  the  magnitude  of  5  is  a  useful  indicator  of  SNR,  which  is  often  used  as  a 
criterion  for  acceptance  or  rejection  of  a  meeisurement. 
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The  relationship  between  members  of  the  set  of  (x,  j/,  u,  v)  distributions  which  result 
in  a  given  {p,q)  distribution  might  be  characterized  as  a  “rotational”  and  a  “magnitude” 
ambiguity.  This  is  because  the  (p,  q)  variables  retain  only  the  phase  and  amplitude  relations 
of  D  to  S.  This  set  contains  many  distributions  that  may  appear  unrelated  to  each  other. 

The  statistical  properties  of  the  MR  have  been  considered  by  Kanter  [1] ,  Tullsson  [2] ,  Seifer 
[3]  and  [4],  and  others.  Most  of  the  previous  studies  have  calculated  the  density  function  of 
real(MR)  and/or  statistical  parameters  of  the  distribution.  Kanter,  in  a  remarkable  study 
[1],  used  sophisticated  techniques  to  obtain  the  density  function  of  p  under  quite  general 
conditions  for  partially- correlated  isotropic  noise,  as  is  considered  here  in  Chapter  5,  but  he 
did  not  discuss  the  joint  distribution  of  p  and  q.  Tullsson  [2]  used  the  method  of  conditional 
probabilities  to  obtain  density  functions  of  both  p  and  q  for  the  same  noise  correlation  model 
under  restrictive  assumptions  regarding  the  mean  values  of  the  sum  and  difference  phasors, 
but  did  not  discuss  the  joint  MR  density.  Both  Kanter  and  Tullson  studied  the  case  where 
the  MR  is  obtained  by  averaging  over  n  observations  and  demonstrated  that  the  variance 
exists  only  forn  >  1.  Seifer  [3]  and  [4]  in  a  remarkable  achievement  has  been  able  to  calculate 
the  mean  value  and  covariance  of  the  MR  under  various  conditions. 

The  joint  distribution  of  p  and  q  is  important  in  systems  that  use  both  the  real  and 
imaginary  parts  of  the  MR  to  determine  target  position.  This  report  considers  the  bivariate 
density  function  of  the  monopulse  ratio  {i.e.,  of  p  and  q)  determined  from  one  trial  n  =  1, 
under  somewhat  more  general  conditions,  and  determines  this  density  function  in  terms  of 
elementary  functions  using  a  straightforward,  elementary  method  to  be  described. 

The  results  of  this  report  are  consistent  with  Kanter’s  and  Tullsson’s  general  conclusions, 
and  in  a  few  easily-checked  special  ceises  their  expressions  are  shown  to  agree  with  those  in 
this  report.  Monte  Carlo  simulations  were  carried  out  to  check  some  of  the  analytically- 
determined  density  functions,  and  to  obtain  the  density  function  in  cases  that  are  difficult 
to  examine  by  analysis.  One  notable  feature  is  the  presence  of  a  “hole”  or  minimum  in  the 
joint  density,  which  occurs  for  some  values  of  the  means  of  the  sum  and  difference  phasors 
and  noise  covariance. 

The  MR  is  a  stochastic  variable  because  of  the  presence  of  various  types  of  noise:  thermal 
or  electronic  noise  in  the  circuitry,  thermal  background  radiation  from  the  sky  which  is 
often  less  important  than  the  other  sources,  diffuse  sea  reflection  noise,  glint,  atmospheric 
variation,  etc.,  all  of  which  contribute  to  the  variance  in  the  distribution  of  the  MR.  Most 
statistical  studies  of  noise,  including  this  one,  assume  that  the  effect  of  the  noise  on  S  and 
D  is  modeled  by  adding  a  Gaussian  error  to  the  true  value.  This  is  a  justifiable  assumption 
for  many  types  of  noise,  while  in  other  cases  it  is  less  appropriate  but  used  anyway  because 
more  detailed  treatments  are  difficult  or  not  feasible.  Noise  also  can  affect  the  mean  values 
S  and  D. 
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Electronic  noise  and  the  diffuse  sea  reflections  add  many  small  random  increments  to 
the  received  signals  and  thus  are  often  taken  to  be  Gaussian.  Atmospheric  variation  and 
glint  are  more  complicated,  and  are  not  representable  as  an  additive  white  Gaussian  noise. 
The  well-known  Swerling  models  specify  the  statistical  distribution  of  the  voltage  amplitude 
factors  for  either  S  or  D,  and  their  correlation  in  time.  The  distribution  of  phases  depends 
critically  on  the  details  of  the  glint  model,  and  while  studies  of  the  effect  of  glint  on  the 
amplitude  of  signals  received  have  been  reported,  to  our  knowledge  very  little  attention  has 
been  given  to  the  effect  on  the  phases. 

Seifer  [3]  has  treated  a  uniformly-distributed  phase  uncertainty,  which  leads  to  the  conclu¬ 
sion  that  whatever  values  S  and  D  would  have  in  the  absence  of  the  phase  uncertainty,  these 
values  become  zero  if  it  is  present.  However,  meaningful  values  of  the  Monopulse  Ratio  D/ S 
are  obtained  because  the  phase  difference  between  S  and  D  is  not  uniformly  distributed  and 
in  fact  peaks  at  a  definite  value.  This  report  describes  cases  in  which  S  and  D  are  allowed 
to  have  arbitrary  values,  which  correspond  to  cases  where  the  phase  uncertainty  does  not 
dominate.  One  apphcation  of  this  more  general  approach  would  be  the  study  of  the  MR 
acquired  from  a  fixed  target. 

The  study  in  this  report  considers  correlated  Gaussian  noise  in  the  sum  and  difference 
channels.  It  is  clear  that  these  assumptions  are  not  valid  in  all  possible  noise  environments. 
While  no  claim  is  made  as  to  the  universal  applicability  of  the  Gaussian  assumption,  the 
results  obtained  on  this  basis  are  useful  in  a  wide  range  of  noise  environments.  The  statistical 
method  used  in  the  study  is  outlined  in  Chapter  2.  Monte  Carlo  methods  of  studying  the 
statistical  distributions  are  discussed  in  Chapter  3.  The  results  for  uncorrelated  isotropic 
noise  are  treated  in  Chapter  4,  while  partially  correlated  isotropic  noise  is  treated  in  Chapter 
5.  The  results  for  arbitrary  Gaussian  noise  are  given  in  Chapter  6.  The  exclusion  of  mea¬ 
surements  when  the  magnitude  of  S  is  below  a  given  threshhold  is  considered  in  Chapter 
7.  An  overall  discussion  of  the  results  is  presented  in  Chapter  8.  Note  that  the  figures  are 
placed  following  the  text  at  the  end  of  each  chapter. 
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CHAPTER  2 

A  CONVENIENT  TRANSFORMATION 


The  development  of  the  MR  density  function  is  facilitated  by  a  well-known  theorem  in 
statistics,  which  states  that  the  density  distribution  of  a  set  of  n  variables  functionally 
related  in  a  one-to-one  manner  with  a  primary  set  of  n  variables  whose  density  function  is 
known  is  obtained  simply  by  substituting  the  functional  relations  into  the  known  density 
distribution  and  multiplying  by  the  Jacobian  of  the  transformation  (see,  for  example,  [5]). 
As  the  density  function  of  only  two  variables,  p  and  q,  as  defined  in  (1.1)  and  (1.2),  is  sought, 
the  immediate  necessity  is  to  define  an  additional  pair  of  dummy  variables  in  order  to  obtain 
a  transformation  that  is  one-to-one  and  onto.  It  is  convenient  to  let  r  =  u  and 

s  =  t;  be  the  dummy  variables.  The  direct  transformation  is  then 


XU  -f-  yv 

P  - 

-f 

while  the  inverse  transformation  is 


g  = 


yu  —  XV 

4.  y2  ’ 


r  =  u. 


S  =  V 


(2.1) 


X  =  pr  ~  qs,  y  =  ps  -{■  qr,  u  =  r,  v  =  s 


(2,2) 


It  is  seen  that  this  transformation  is  indeed  one-to-one  over  the  infinite  four-dimensional 
space.  The  Jacobian  of  the  transformation  is  given  by 


d{x,y,u,v) 

r 

—s 

P 

-9 

s 

r 

q 

p 

d{p,q,r,s) 

0 

0 

1 

0 

0 

0 

0 

1 

2  ,  2 
=  r‘  s 


(2.3) 


This  formulation  is  convenient  for  the  case  where  there  is  a  threshhold  Rq  on  the  sum  phasor 
modulus.  If  the  joint  distribution  of  u  and  v  is  modified  to  exclude  values  <  Rq 

(which  would  require  renormahzation  of  the  density  function),  the  corresponding  marginal 
distribution  of  p  and  q  is  obtained  by  restricting  r  and  s  to  values  -f-  >  Rl 

Seven  cases  will  be  considered  in  order  of  increasing  complexity.  Chapter  4  considers  the 
first  four  cases  in  which  the  noise  in  the  sum  and  difference  channels  is  assumed  uncorrelated 
and  isotropic.  In  case  I  the  means  of  the  sum  and  difference  phasors  are  taken  to  be  zero. 
In  case  II  the  same  assumptions  are  made,  except  that  the  sum  channel  is  allowed  to  have 
a  specified  nonzero  mean.  Case  III  is  the  same  as  Case  II,  except  that  only  the  difference 
phasor  is  allowed  to  have  a  specified  nonzero  mean.  Case  IV  is  the  same  except  that  the 
phasors  of  both  channels  are  allowed  to  have  specified  nonzero  means. 
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Tullsson  [2]  has  pointed  out  that  it  may  not  be  reasonable  to  neglect  the  noise  correlation 
between  the  sum  and  difference  channels  and  considered  a  correlation  between  the  real  parts 
and  between  the  imaginary  parts  of  the  two  channels.  Case  V  treats  the  simple  case  where 
all  the  mean  values  are  zero,  while  Case  VI  considers  general  means.  Finally,  Case  VII 
considers  an  arbitrary  fourth-order  covariance  matrix  as  well  as  arbitrary  means. 

It  is  well  known  that  the  density  function  of  the  monopulse  ratio  has  a  mean  value  but 
infinite  variance.  That  is,  the  integral  defining  the  variance  diverges  because  of  the  long  tails 
of  the  density  function.  Of  course,  in  any  practical  situation  the  monopulse  ratio  would  be 
bounded,  and  the  variance  could  be  computed.  The  essential  consideration  is  the  manner 
in  which  the  density  function  should  be  modified  to  incorporate  such  limitations  on  the 
permissible  values  of  the  monopulse  ratio  in  order  to  arrive  at  a  satisfactory  value  for  the 
variance.  Tullsson  [2]  suggests  a  “threshholding”  idea,  according  to  which  measurements 
having  a  sum-channel  magnitude  smaller  than  a  given  quantity  are  discarded.  This  is  equiv¬ 
alent  to  considering  a  modified  Gaussian  density  where  a  circular  region  in  the  (r,  s)  plane 
is  excluded.  The  resulting  monopulse- ratio  densities  are  difficult  to  calculate  in  general,  and 
are  presented  only  in  a  few  simple  cases. 

Let  x,y,  M,  u  be  the  mean  values  of  x,y,u,v,  and  let  (Txy,<Txu--  be  the  components  of  the 
covariance  matrix.  “Isotropy”,  as  in  Cases  I  through  VI  in  Chapters  4  and  5,  is  understood 
to  mean  that  there  are  the  relations 


defining  the  quantities 
is  understood  to  mean 


(Xxx  —  (^yy  — 


(^uu  —  (^vv  — 


(2.4) 


ax  and  or„.  “Uncorrelated”,  as  in  Cases  1  through  IV  of  Chapter  4, 
that 


(^xy  —  ^xu  —  ^xv  —  (^yu  —  (^yv  —  ^uv  —  0 


(2.5) 
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CHAPTER  3 

NUMERICAL  GENERATION  OF  RANDOM  VECTORS 

It  is  useful  to  generate  series  of  random  vectors  or  linear  arrays  having  a  specified  joint 
probability  density.  The  density  function  of  interest  here  is  the  joint  Gaussian  distribution 
of  the  four  variables,  x,  y,  u,  v.  However,  the  method  to  be  discussed  works  for  any  number  n 
of  dimensions.  First,  obtain  n  independent  normal  variates,  o,-,~  N{ai,l),i  —  l,...,n,  and 
construct  the  desired  n-variate  by  a  linear  transformation 

X=^QA  (3.1) 

where  X  —  (xi,X2,  ...,x„)^  is  the  desired  n-variate,  A  =  (Q;i,a2,  ...,a„)^,  and  Q  is  the  n  by 
n  transformation  matrix  to  be  determined.  The  a  variables  thus  have  specified  means,  unit 
variance,  and  are  independently  Gaussian  distributed.  The  desired  covariance  matrix  of  X 
is 

S  =  E[{X  -  X){X  -  Xf]  =  E[XX'^]  -  XX'^  =  QE[AA^]Q^  -  XX'^  (3.2) 

where  the  overbar  indicates  the  mean  or  expected  value.  According  to  the  stated  statistical 
properties  of  the  a  variables,  it  is  seen  that  •. 

E[AA^]  =  AA^  +  I  (3.3) 

where  /  is  the  n  by  n  identity  matrix.  Equation  (3.2)  then  becomes 

E  =  Q{AA^  -h  I)Q'^  -  XX^  =  QAA^Q'^  +  QQ'^  -  XX'^ 

=  XX^  +  QQ^  -  XX'^  =  QQ^ 

Since  S  is  symmetric,  a  symmetric  solution  Q  of  (3.4)  usually  can  be  obtained  by  just  taking 
the  square  root  of  the  matrix  E.  There  is  usually  a  multiplicity  of  solutions,  from  which  any 
one  real  square  root  matrix  will  provide  the  desired  result.  The  means  of  the  a  variables  are 
given  by 

A  =  Q-^X  (3.5) 

Some  numerical  computing  packages  (e.g.,  Matlab)  have  a  matrix  square  root  routine, 
which  in  many  well-behaved  cases  determines  a  real  square  root  of  a  real  square  matrix.  If 
such  a  routine  is  not  available,  an  eigenvalue  routine  can  be  used  as  follows.  While  there 
may  be  more  elegant  or  more  general  procedures,  the  following  meets  the  immediate  need 
and  is  simple.  Let  A  and  v  be  an  eigenvalue  and  eigenvector  of  Q.  Then, 

Qv  =  Xv 
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Multiplying  on  the  left  by  Q  gives 

QQv  =  Tiv  =  XQv  =  X^v  (3-7) 


showing  that  the  eigenvalues  of  Q  are  just  the  square  roots  of  the  eigenvalues  of  S,  and  their 
eigenvectors  are  the  same.  Let  V  be  the  n  by  n  matrix  whose  columns  are  the  eigen vectois 
of  E,  and  let  be  the  diagonal  matrix  whose  elements  are  the  eigenvalues  of  S.  The 
matrix  L  is  diagonal  and  has  elements  A  which  are  the  square  roots  of  the  elements  of  L  . 
Then 

(3.8) 

QV  =  VL 

from  which 

Q  =  VLV-^  (3.9) 


is  the  desired  solution. 


These  results  are  used  to  generate  Monte  Carlo  stochastic  samples  of  the  (x,  y,  u,  v)  space 
having  the  desired  density  function,  from  which  the  resulting  monopulse-ratio  quantities 
{p,q)  are  obtained  from  the  relations  (3.2)  and  plotted  as  a  scatter  diagram.  Such  diagrams 
provide  a  useful  demonstration  of  the  validity  of  the  density  expressions  obtained  by  analysis. 
The  method  is  also  useful  for  empirical  determination  of  density  functions  in  cases  where 
analytic  solutions  are  not  available.  It  will  be  noted  in  subsequent  chapters  that  the  marginal 
density  function  for  p  alone  (or  q  alone)  is  often  difficult  to  obtain  analytically,  in  which  case 
the  Monte  Carlo  method  can  provide  a  useful  estimate.  However,  analytic  solutions  are 
preferable  because  their  dependence  on  the  various  parameters  is  immediately  visible. 
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CHAPTER  4 


UNCORRELATED  ISOTROPIC  NOISE 


The  simplest  although  not  the  most  realistic  case  is  where  the  noise  covariance  matrix  is 
diagonal.  In  accordance  with  the  assumed  isotropy  the  covariance  noise  matrix  is  assumed 
to  be 


( <Pi 

0 

0 

0\ 

0 

0 

0 

0 

0 

0 

Vo 

0 

0 

^u/ 

where  <Ty  is  taken  equal  to  ax,  and  a^  equal  to  cr„.  Four  subcases  are  considered  by  first 
taking  the  means  of  all  the  primary  variables  (a;,  y,  u,  v)  to  be  zero  amd  finally  allowing  these 
means  to  have  arbitrary  values. 


CASE  I  -  Zero  Means 


The  means  x,  y,  u,  v  are  taken  to  be  zero.  The  density  distribution  of  the  variables  x,  y,  u,  v 
is  accordingly  given  by 


1 


Pr,y,uA^,y,U,v)^j-^^exp 


x^  +  y^  -f 

2cr2 


2al 


(4.2) 


According  to  the  statistical  theorem  stated  in  Chapter  2,  the  density  function  of  the  variables 
p,  q,  r,  s  is  given  by 


Pp,q,T,s{P^  9)  «)  =  Px,y,u,v{p'^  “  9^7 +  qr,  r,  s) 


d{x,y,u,v) 


47r2o-2cr2 


47r2(T2(72 


d(p,y,r,s) 
(pr  —  qs)^  +  (ps  +  qr)^ 


2al 


2al 


(4.3) 


exp|  —  (r^  +  s^) 


2ai 


2al 


The  marginal  density  of  the  p,  q  variables  is  found  by  integrating  the  above  expression  over 
the  entire  (r,  s)  plane.  The  substitution 


r  =  R  cos  9 
s  =  Rsin  9 


(4.4) 


gives 


1  /•23r  r  / 

exp[-fi^( 


,2  f  pI±^  J_ 

2^^  2a^, 


d9 


(4.5) 
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The  result  is 


(4.6) 


where  7  is  written  for  the  ratio  (Ti/cTu  of  the  standard  deviations  (see  Appendix  for  details). 
This  result  has  been  presented  by  Tullsson  [2]  under  somewhat  more  general  conditions.  It 
is  noteworthy  that  the  above  joint  density  distribution  depends  not  on  the  two  individual 
variances  but  only  on  their  ratio. 

The  density  function  (4.6)  can  be  expressed  in  more  compact  form  containing  no  param¬ 
eter  by  defining 

p'  =  ph,  ^p,g  =  7^^P,9  (4-7) 

Then  .. 

Pp,gip'^  Q')  =  ^(^,2+^/2  + 1)2 

Essentially  the  parameter  7  has  disappeared  through  appropriate  scaling.  This  function  is 
shown  graphically  in  Figures  4-1  and  4-2  at  the  end  of  this  chapter. 


It  is  easy  to  find  the  marginal  density  of  p  by  integrating  (4.6)  with  respect  to  q.  We  have 

>^2  dq 

U^212 

(4.9) 


[p2  4.  ^2  ^2]2 


A  similar  compact  form  of  the  marginal  density  (4.9)  is  obtained  through  the  scaling  (4.7) 
along  with  the  definition 

Pi  =  lPp  (4.10) 


giving 


p;ip')  = 


1 


(4.11) 


7r(p'2  + 1)3/2 

This  function  is  shown  in  Figure  4-3.  There  are  similar  expressions  for  the  marginal  density 
of  q.  It  is  readily  demonstrated  that  the  integral  of  Pp{p)  in  (4.9)  over  the  entire  p  axis  is 
unity.  The  mean  value,  p,  is  zero  because  (4.9)  is  an  even  function  of  p  while  the  integral 
defining  p  converges.  However,  the  variance  of  p  is  indeterminate  as  it  is  expressed  by  a 
divergent  integral. 


CASE  II  -  Zero  Mean  of  the  Difference  Channel  Phaser 


The  covariance  matrix  (4.1)  and  means  x  and  y  of  zero  are  assumed,  while  u  and  v  can 
be  assigned.  The  density  function  for  the  primary  variables  is 


Px,y,u,vix,y,U,v) 


x^  +  y"^  {u  —  u)"^  +  {v  —  vY 


(4.12) 
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The  density  function  of  the  secondary  variables  is  then  given  by 


(pj —  qs)^  +  {jps  +  qrY  {r  —  u)^  +  (s  —  u) 


21 


47rW2 


-»-( 


2ctI 

p2  4.  g2  ^  1 


2(7.2 


where 


U  = 


2(72  2(7.2 


2(7.2 


(r^  +  3")  + 


u 

ur  4-  vs 


(4.13) 


(4.14) 


To  find  the  marginal  distribution  of  (p,q)  the  substitution  (4.4)  is  again  used  to  obtain 

roo  r  /p2^g2  I 


e-u  foo 

PpAp^Q)  =  I  P^dRexp  - 


2(72  +  2(7,2 


/  exp 
Jo 


[  R{ucos6  +  usin^) 


L 


(4.15) 


de 


The  details  of  the  integrations  are  given  in  the  Appendix.  The  inner  integral  is 


l‘2jt 

/  exp 
Jo 


R(ucos9  +  usin^) 


d6  =  277 /o 


S'* 


u2  +  u2 


(4,16) 


where  Iq  is  the  Modified  Bessel  Function  of  zero  order.  The  outer  integral  is  then 


Pp,q{P^q)  = 


,-U  roo 


2Tr  a-‘2  a. 


d 


R^dR  exp 


-( 


P  +g  ,  J_ 
2(72  2(7.2 


R^ 


h  ^^\/u2  +  u2^  (4.17) 


which  becomes 


Pp,q{P^(l)  = 


p^(l  +p)e^ 

777/272 


where  p  stands  for  the  quantity 


p  = 


t/72 


U 


p2  -I-  g2  4.  ^2  pi2  4.  g(2  4.  1 


(4.18) 


(4.19) 


It  is  reasily  demonstrated  (see  Appendix)  that  the  total  volume  under  the  joint  density 
function  (4.18)  is  unity,  and  it  is  seen  that  (4.18)  reduces  to  (4.6)  of  Case  I  if  77  —>  0.  The 
compact  form  of  (4.18)  is 

(4.20) 

with  the  definitions  (4.7)  and  (4.19).  This  function  is  shown  in  Figure  4-4  for  various  values 
of  U. 


The  marginal  density  of  p  is 


p-t/  roo 

J  P^{^+p)e''dq  (4.21) 
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As  pi  depends  on  q  it  is  convenient  to  use  pi  as  the  variable  of  integration.  The  details  aie 
given  in  the  Appendix.  After  the  transformation  the  density  function  is 


p{\  +  pL)e>^dpi 


which  finally  leads  to 


Pp{p)  = 


E 


V'I/72/x  -  (p2  +  72)/x2 


(2j  +  l)!  /  fV 


or  to 


2(p2  +  72)3/2  ^  [;!]3  \^4(p2  +  ^2)y 

PVVl  -  V  ^  V 

)- 2(p«  +  1)3/2  ^  \j\]3  V4(p^  +  1)/ 


In  terms  of  standard  functions  these  expressions  are  equivalent  to 


Ppip)  = 


r 


2(p2  + 72)3/2 

t/72 


or  to 


P^ip')  = 


p2  4.  ^2 

1 

2(p'2  +  1)3/2 


exp 

/o 


_Uf2£+J^\ 
2  Vp^+7^  /. 
U'r‘^ 


X 


[/  _u±_\  J  _  Jlf_T  (  Y 

A  '^P^  +  7^/  °\  2(p2+72)y  ^24.72  ^  2(p2+72)y 

u  /2p'2  +  r 


exp 

h 


2  V  P'2  +  1 

U 


2(p'2  +  l)) 


U 


-h 


V 


(4.22) 


(4.23) 


(4.24) 


(4.25) 


(4.26) 


p'2  +  i;‘"v  2(p'2  +  i);  p'2  +  r*v  2(p'2  +  i);j 

where  /q  and  h  are  Modified  Bessel  functions.  This  density  function  is  shown  in  Figure  4-5 
for  various  values  of  U. 

Again  it  is  noted  that  p  =  0  and  the  variance  of  p  is  infinite.  The  expressions  (4.25)  and 
(4.26)  reduce  to  (4.9)  of  Case  I  when  t/  0  because  /o(0)  =  1  and  /i(0)  =  0. 


CASE  III  -  Zero  Mean  of  the  Sum  Channel  Phaser 

The  covariance  matrix  (4.1)  and  means  u  and  v  of  zero  are  assumed,  while  x  and  y  can 
be  assigned.  For  this  case  the  density  function  for  the  primary  variables  is 


(x  -  a;)2  +  {y-  y)^  _ 

2cr2  2(t2 


Fi,y,M,t,(x,y,u,u)  —  ^^2(j2^2 

The  density  function  of  the  secondary  variables  is  then  given  by 

{pr  -qs-  x)2  -j-  [ps  +  qr  -  p)2  _  r2  +  ^2 


(4.27) 


r2  +  52 


/p,5,r,s(p,9,r,s)  -  ^^2fj2cr2 


x''u 


2(71 


2(71 


r2  -1-  s2 
47r2(T2cr2 


exp 


-X 


(r2  +  32)(p2  4.  g2)  _  2p(rx  +  sy)  +  2q{sx  -  ry)  _  +  s 


21 


2^2 


2<^2  j 


(4.28) 
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where  X  stands  for  the  quantity 


X  = 


2(jI 


(4.29) 


The  marginal  distribution  of  p  and  q  is  again  found  by  means  of  the  substitution  (4.4),  to 
obtain 


-A  foo 


2^2  2<t2 


X 


(4.30) 


I  dO  exp  —  {\px  4-  qy]cos6  +  [py  —  yx]sin0) 

Jo 

The  inner  integral  is  evaluated  in  same  way  (4.17)  was  obtained  from  (4.15).  The  result  is 


I  +  2^)]  'o  +  + 

(4.31) 

This  integration  is  performed  in  the  same  way  that  (4.18)  was  obtained  from  (4.17),  giving 


Pp,qiP^^) 


7^p^(l  +  p)e^  ^ 

7rA’2(p2  q.  q2y 


(4.32) 


or 


where  p  is  now  defined 


p^(l  +  p)eP  ^ 

7rX2(p/2  4.  ^/2)2 


_  (p^  +  q^)X  _  jp'^  +  q’^)X 

^  p2  q.  ^2  4.  ^2  P'2  q.  qi2  _|.  j 


(4.33) 


(4.34) 


Again,  p  =  0  because  (4.32)  is  an  even  function  of  p  and  the  integral  defining  p  converges. 
The  integral  defining  the  variance  does  not  converge,  so  again  no  finite  variance  exists.  The 
expression  (4.32)  reduces  to  (4.6)  of  Case  I  when  A”  — ^  0.  The  marginal  distribution  of  p  is 
somewhat  more  cumbersome  to  obtain  than  in  the  previous  cases,  and  has  not  been  devel¬ 
oped.  Figure  4-6  shows  this  joint  density  function  for  various  values  of  X.  It  is  noted  that 
there  is  a  minimum,  or  depression,  at  the  origin  for  values  of  X  greater  than  approximately 
unity. 


CASE  IV  -  General  Means 


For  this  case  the  covariance  matrix  (4.1)  is  assumed  while  the  means  of  all  four  primary 
variables  x,y,u,v  are  specified  to  be  x,  y,  u,  v.  The  density  function  of  the  primary  variables 
is  accordingly  given  by 


P t,y,u,v{x,y,  U,  v)  — 


47r2(r2(72 


exp 


(x  -  x)2  +  {y  -  yr)  {u  -  uy  +  {v  -v) 


,-.^2^ 


(4.35) 
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The  corresponding  density  function  for  the  secondary  variables  is 


Pv,q,T,s{P,q,'r,s)  = 


'x'^u 


exp 


(pr  -qs-  xf  +  {ps  +  gr  -  yf  _  (r  -  +  {s  -  v)'^ 

2(tI 

fp^  + 


2^^ 


^^yr^W)-X-U 


X 


'xp  +  yq  U 

<%  I  A 


,  fyp-xq  ^  V 

r  +  - - 1- 


(4.36) 

The  integration  can  be  carried  out  in  the  manner  that  was  used  to  obtain  (4.32)  from  (4.28). 

The  result  is  ^ 

,  -yHl  +  n)  el-X-" 

,2  +  .,2)2 


or 


p,  (I  ,,  ^  (l+)l)e<‘  ^  ^ 
’“I  I  j(p/2+,/2  +  ))2 


(4.37) 

(4.38) 


where  here  p  is  redefined  to  be 


X{p‘^  +  q^)  +  'yjWp  +  Zq)  +  7^(7  _  X{p'‘^  +  q'^)  +  Wp'  4-  Zq'  +  U 


p2  j^q2  ^  ^2 


P'2  +  q>^  +  l 


with  the  parameters 


W  = 


ux  +  vy 


Z  = 


uy  —  vx 


(4.39) 


(4.40) 


(TxO'u  (XxO'xi 

which  depend  only  on  the  density  function  of  the  primary  variables,  being  defined  for  conve¬ 
nience.  When  X  =  p  =  0,  (4.37)  reduces  to  (4.18)  of  case  II.  When  u  =  u  =  0,  (4.37)  reduces 
to  (4.32)  of  Case  III.  The  density  (4.37)  is  not  an  even  function  because  of  the  presence  of  p 
and  q  to  the  first  power  in  the  expression  (4.39)  for  p.  Consequently  p  and  q  do  not  vanish 
as  they  did  in  previous  cases.  Figures  4-7  and  4-8  show  the  function  (4.37)  for  the  case 


x  =  2,  y  =  2,  u  =  5,  v  =  0,  cr^  =  1,  cr„  =  20 


(4.41) 


Figure  4.9  shows  various  perspectives  of  the  same  function.  A  Monte  Carlo  simulation  was 
carried  out  by  selecting  random  values  of  the  primary  variables  from  the  density  function 
defined  by  (4.35)  and  calculating  the  corresponding  complex  monopulse  ratios.  These  are 
plotted  on  the  complex  plane  as  a  scatter  diagram  and  shown  in  Figure  4-10  as  an  experimen¬ 
tal  validation  of  the  derived  joint  density  function  (4.37).  Thus,  Figure  4-9  is  an  empirical 
counterpart  to  the  joint  density  function  depicted  in  Figures  4-7,  4-8,  and  4-9. 
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FIGURE  4-1.  JOINT  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  MEANS 


FIGURE  4-2.  JOINT  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  MEANS 


P(p,q)»qamma~2  P(p)»qamma 


FIGURE  4-3.  REAL  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  MEANS 


sqrt(p-2+q-2)/gomma 

FIGURE  4-4.  JOINT  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  DIFFERENCE  MEAN 
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FIGURE  4-5.  REAL  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  DIFFERENCE  MEAN 


Wonopulse  Density  for  Cose  ill,  X  —  0,1 . 8 


FIGURE  4-6.  JOINT  MR  DENSITY,  UNCORRELATED  NOISE,  ZERO  SUM  MEAN 
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FIGURE  4-9.  JOINT  MR  DENSITY,  UNCORRELATED  NOISE,  GENERAL  MEANS 
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CHAPTER  5 

PARTIALLY  CORRELATED  ISOTROPIC  NOISE 


Noise  effects  in  the  sum  and  difference  phasors  are  correlated  because  they  are  linear 
combinations  of  the  same  antenna  lobe  outputs.  We  shall  consider  the  noise  covariance 
matrix  suggested  by  Tullsson  [2],  according  to  whom  the  assumption  of  real  antenna  patterns 
gives  a  correlation  between  the  real  parts  and  between  the  imaginary  parts  but  not  crosswise. 
Accordingly,  the  covariance  matrix  can  be  represented 


/ 

0 

p(Tx  O'-ii 

0  > 

( 

7^ 

0 

Pl 

0  \ 

0 

0 

pCTx^u 

=  (t2 

0 

A 

0 

Pl 

p(T  xCT  u 

0 

0 

^  U 

P7 

0 

1 

? 

0 

p(Tx 

0 

/ 

\ 

0 

Pl 

0 

1  / 

S  = 


where  p  is  a  correlation  coefficient.  The  determinant  and  inverse  of  (5.1)  are 


1^1  =  -  P^) 


2\2 


= 


1 


1 

0 

-Pl 

0 


0 

1 

0 

-PI 


-PI 

0 

0 


0 

-pl 

0 

7^  / 


(5.1) 


(5.2) 


CASE  V  -  Zero  Means 


When  the  means  x,  y,  u,  v  are  taken  to  be  zero,  the  density  function  of  the  primary  vari¬ 
ables  is 

1 


Px,y,u,v{^,y,u,v)  = 


X 


exp 


47rV2<T2(l  -  p2) 

-  2pax(TuXU  -f  crlu^  +  aly'^  -  2pax(Tuyv  -|-  cr2t;2 


The  corresponding  density  function  of  the  secondary  variables  is  then 


(5.3) 


Pp,?,r,s(p,g,r,s)  = 


r2  -|-  ^2 


•X 


47rV2(T2(l  -  p2) 

-  g-s)^  -  pcTif^uripr  -  qs)  +  alr^  -f  cr^(p.s  -h  qr)^  -  2p(TxCrus{ps  -f  qr)  -f 


exp 

After  some  rearrangement  this  becomes 


-  P^) 


) 


Pp,q.rAP^  9^  A  = 


47rV2o-2(l  -  p2) 


exp 


(p^  +  q'^  -  2p7p  +  7^)(r2  -|-  g2) 
2tr2(l  -  p2) 


(5. 
(5.5) 
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The  marginal  distribution  of  the  variables  {p,q)  is  obtained  as  in  (4.6)  with  the  result 

P  (  \  7^(1 /Kg) 

This  function  reduces  to  the  joint  density  (4.6)  if  p  is  set  to  zero.  Equation  (5.6)  is  illustrated 
in  Figures  5-1  and  5-2  for  various  values  of  p.  The  compact  form  is 


p!  f  I  l\ _ ^  _ 

9  +  gl2  _  2pp'  +  1)2 


The  marginal  distributions  are 


/OO 

^P,?(PW)<^9  = 

-00 


7^(1  - 


and 


2(p2  —  2p'yp  +  72^3/2 

7^(1  -p2) 


PM  =  P..M<l)dp  = 

The  corresponding  compact  forms  are 

pi(n')  =  (1-^^) 

2(p'2_2pp'-M)3/2 


and 


p.(V)  = _ (w!) _ 

>  2[?'2  +  (1  _  ^2)]3/2 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


(5.11) 


The  density  functions  (5.10)  and  (5.11)  are  illustrated  in  Figures  5-3  and  5-4  for  various 
values  of  p.  It  is  seen  that  q  =  0.  If  the  substitution  p  =  t  +  'jp  is  made  in  (5.8),  it  is  seen 
that  f  =  0,  which  shows  that  p  =  p^.  The  variances  of  both  p  and  q  diverge  as  in  the  case 
of  uncorrelated  noise. 


CASE  VI  -  General  Means 


When  the  means  of  the  primary  variables  can  be  specified,  the  joint  density  function  of 
the  secondary  variables  is 


r2  -f  52 


-exp 


A(r2  -|-  s2)  -f  Br  +  Cs  +  D 

2(t2(1  -p2) 


47r2(T2o-2(l  -  p2) 

where  A,  C,  and  D,  and  their  primed  counterparts  are  defined  by 

j^A!  =  A  =  {p^  +  q^)  -  2pjp  -}-  7^ 
ja^B'  =  B  =  -2{xp  +  yq)  +  2p'y{x  +  up  +  vq)  -  2'^'^u 
'jaxC  =  C  =  2{xq  -  yp)  -  2pj{uq  -  y  -  vp)  -  2'^‘^v 
alD'  =  D  =  (x"^  +  y^)  -  2p'y{vy  +  ux)  -t-  7^(fi^  -f  v"^) 


(5.12) 


(5.13) 
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all  being  independent  of  r  and  s.  A  feature  of  (5.12)  which  facilitates  further  development 
is  the  absence  of  a  product  rs  either  inside  or  outside  the  exponential  factor.  To  get  the 
marginal  bivariate  density  of  p  and  q,  the  integration  with  respect  to  r  and  s  is  again  carried 
out  with  the  aid  of  the  transformation  (4.4)  to  obtain 


Pp,g{P^<l) 


47r2<T2(T2(l  -p2) 


f 


R^exp 


AR?  +  D 


L  2^2(1 -p2) 


dR  X 


/  exp 
Jo 


R{Bcos$  +  Csin^) 


2<t2(1  -  P‘2) 


(5.14) 


de 


The  inner  integration  can  be  carried  out  in  same  way  that  (4.19)  was  obtained  from  (4.15). 
The  result  is 


Pp,q{P^(l) 


27r(r2cr2(l 


2(1  -  p‘‘)  I 


/*oo 

r  AR^  +  D  ] 

lo 

RVB^  +  (72 

/  i?'*exp 
Jo 

.  2o-2(l-p2)_ 

2(t2(1  -p2)^ 

dR  (5.15) 


The  last  integration  can  be  carried  out  in  the  same  way  that  (4.18)  was  obtained  from  (4.17), 
giving 


PpAp^q) 


7^(1  -P^) 

■kA2 


B^  +  C^  ■ 
8X<t|(1  -  ,,2), 


exp 


/  4AD-B^-C^\ 
V  8At72(l-p2)  J 


(5.16) 


It  can  be  shown  that  (5.16)  reduces  to  (4.37)  if  p  — +  0.  An  example  of  this  joint  density  func¬ 
tion  is  shown  in  Figures  5-5  and  5-6  for  the  values  (4.41)  and  p  =  .99.  Various  perspectives 
of  this  function  are  shown  in  Figure  5-7.  Comparing  Figures  5-5,  5-6  and  5-7  for  p  =  .99 
with  Figures  4-7,  4-8,  and  4-9,  for  p  =  0  reveals  the  elFect  of  the  correlation  between  the 
sum  and  difference  phasors  on  the  resulting  joint  MR  density.  In  particular,  it  is  noted  that 
the  correlated  case  has  a  deeper  and  flatter  ‘hole’  or  minimum.  A  Monte  Carlo  simulation 
was  carried  out  in  the  manner  of  Figure  4-10  and  the  results  are  presented  in  Figure  5-8. 
Comparing  the  two  scatter  diagrams  shows  empirical  confirmation  that  the  correlation  leads 
to  a  deeper  minimum  in  the  joint  density. 


A  compact  form  of  (5.16)  is  obtained  by  defining 

x'  =  i/o-i,  y'  =  yfax,  u'  =  u/(7u,  v'  =  V  j 


With  these  definitions  the  primed  quantities  defined  in  (5.13)  become 

A'  =  p'2  +  q>2  _  2pp'  +  1 

B'  =  —2{x'p'  -I-  y'q')  4-  2p{x'  -f  up'  -t-  v'q')  —  2u' 
C'  =  2{x'q'  —  y'p')  —  2p{u'q'  —  y'  —  v'p')  —  2v' 
D'  =  x'^  +  y'^  —  2p{v'y'  -t-  u'x')  +  u'^  -i-  v'^ 
These,  along  with  (4.7)  allow  (5.16)  to  be  written  in  the  form 


Pi,,g{P^9) 


IzZ. 

7rA'2 


B'2-pC'^  ■ 
8A'(l-p2) 


exp 


{  SA'{l-p^)  ) 


(5.17) 


(5.18) 


(5.19) 
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Cose  V - rho  =  .9 


FIGURE  5-1.  JOINT  MR  DENSITY,  PARTIALLY  CORRELATED  NOISE 


p/gommo 


FIGURE  5-2.  JOINT  MR  DENSITY  FOR  VARIOUS  VALUES  OF  THE  CORRELATION  COEFFICIENT 
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FIGURE  5-3.  REAL  MR  DENSITY  FOR  ZERO  MEANS,  FOR  VARIOUS  VALUES 
OF  THE  CORRELATION  COEFFICIENT 


q/ gommo 


FIGURE  5-4.  IMAGINARY  MR  DENSITY  FOR  ZERO  MEANS,  FOR  VARIOUS 
VALUES  OF  THE  CORRELATION  COEFFICIENT 
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5-7 


FIGURE  5-7.  JOINT  MR  DENSITY,  PARTIALLY  CORRELATED  NOISE,  GENERAL  MEANS 
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FIGURE  5-8.  MONTE  CARLO  SCATTER  DIAGRAM  OF  MR  FOR  PARTIALLY  CORRELATED  NOISE 
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CHAPTER  6 

GENERAL  CORRELATED  NOISE 


CASE  VII  -  General  Means 


Here  the  basic  assumption  of  gaussian  density  for  the  primary  variables  x,  y,  u,  v  is  re¬ 
tained,  but  the  covariance  and  mean  values  are  arbitrary.  The  general  four-variate  gaussian 
density  function  is 

Px,j,.a.«(a;,i/,u,n)  =  j  exp  -  ^)]  (6-1) 

47r'^  |Sj  ' 


where  ^  stands  for  the  column  vector  containing  the  primary  variables  x,y,u,v.  Let  the 
inverse  of  the  covariance  matrix  be  expressed  by 


^  ai\ 

a\2 

013 

014  \ 

a\2 

022 

023 

024 

ai3 

023 

033 

034 

\«14 

024 

034 

044  / 

(The  covariance  matrix  and  its  inverse  are  taken  symmetrical  without  loss  of  generality.) 


The  transformation  (2.2)  is  linear  in  r  and  s,  and  can  be  written  as 


^  =  ar  -h 


(6.3) 


where  a  and  /?  represent  vectors 


( 

f-9\ 

P 

0 

VO/ 

V  1/ 

(6.4) 


which  are  independent  of  r  and  s.  With  the  help  of  these  vectors  the  exponent  in  (6.1)  can 
be  written 


—  ^a)r^  —  (a^S  ^/3)rs  —  ^(/S^S  ^ 

+  (a^S-‘|)r  +  (/?’'S-"|)s  - 


=  —{Ar^  -}-  Brs  4-  Cs^  +  Dr  Es  +  F) 
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thereby  defining  the  quantities 

0  0 

A  =  ^anp  +  ai2Pg  +  ai3P  +  ^0229  +  «239  +  j«33 

B  =  —a\\pq  +  a\2{p^  -  q^)  -  a^q  +  aup  +  a22Pq  +  a2ZP  +  «249  +  o-U 

C  =  \anq^  -  a\2pq  -  auq  +  \a22P^  +  a2AP  +  §044 

D  =  -auxp  -  a\2{yp  +  xq)  -  an{up  +  x)-  auvp  -  a22yq  -  a,2z{uq  +  y) 

—  a2Avq  —  a33U  —  (6.6) 

E  =  auxq  -  ai2{xp  -  yq)  +  aisu^  —  ai4{x  —  vq)  —  a22yp  —  a23«P 

-  a24{vp  +  y)  —  a34«  —  duv 

F  =  \aux‘^  +  a\2xy  +  ais^u  +  auxv  +  ^022^^  +  a23yu  +  ^24^ 

+  |a33fi^  +  a34fih  +  ia44h^ 

which  are  independent  of  r  and  s.  The  joint  marginal  density  function  of  p  and  q  will  now 
be  developed.  Since  it  is  not  convenient  to  integrate  (6.1)  by  means  of  the  substitution  (4.4), 
the  integration  will  be  carried  out  directly  in  the  variables  r  and  s. 

The  marginal  density  function  of  the  variables  p,9,r  is 

Pp,g,TiP,q,r)  =  — ;; - T7X  /  {r^  +  s^)  exp[-Cs^  -  {E  +  Br)s  -  {Ar'^  +  Dr  +  F)]ds  (6.7) 

47r2  J-00 


which  becomes 

Pp,g,T{p,  Qi  ^) 
where  the  substitutions 


(7rC)-3/2e-<^ 

svisi 


[l+2C{r^  +  r^)] 


E  +  Br 
2C 


G=-Ar^  +  Dr  +  F-  Cr^ 


have  been  made  for  convenience. 


(6.8) 


(6.9) 


Equation  (6.8)  is  meaningful  only  if  (7  >  0,  the  validity  of  which  is  evident  from  (6.5), 
according  to  which 

C  =  (6.10) 

As  S  (and  consequently  also  S~^)  is  positive  definite  in  any  significant  situation,  (6.10)  must 
be  positive  for  any  nonvanishing  /5,  and  in  particular,  for  the  /?  given  in  (6.4). 


The  shape  of  the  three-dimensional  density  function  (6.8)  can  be  counterintuitive.  Con¬ 
sider  the  case  where 


/ 

3.0 

-0.5 

1.0 

0.5  \ 

2\ 

- 

-0.5 

1.0 

5.0 

-0.5 

-0.5 

4.0 

1.0 

0.0 

9 

V 

0.5 

1.0 

0.0 

5.0  y 

^  ly 

(6.11) 
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It  is  difficult  to  display  this  three-dimensional  function  by  means  other  than  presenting 
two-dimensional  sections  accross  it.  Figures  6-1  through  6-4  illustrate  a  few  of  such  sections. 


The  integration  of  (6.8)  over  r,  to  obtain  the  marginal  density  function  of  p  and  q,  gives 


Pp,qiP^<l) 


Ke-^ 

27rA2yri^ 


(6.12) 


where 

A  =  52  -  4AC 

M  =  F  +  A-\AE^-BDE  +  CD^) 

K  =  -f  B^E^  -  ABODE  -1-  40^ -  AABDE  +  AA^E‘^  -f  %A^C 

-  2B^C  -}-  SAC^  -  2AB^ 


The  steps  are  outlined  in  the  Appendix.  It  is  of  interest  to  consider  the  form  of  the  depen¬ 
dence  of  Pp^q  on  p  and  q.  The  quantities  A  and  K  are  polynomials  of  order  four  and  six  in 
{p,q),  while  M  is  a  rational  expression  of  order  four  over  four. 


Equation  (6.12)  is  meaningful  only  if  A  <  0.  The  following  discussion  does  not  pretend 
to  be  a  rigorous  proof  that  this  is  always  the  case,  but  tries  to  make  it  seem  reasonable.  In 
the  same  way  that  C  was  expressed  as  the  quadratic  form  (6.10),  A  and  B  can  be  similarly 
expressed  by 

B  =  (6.14) 

The  quantity  A  is  now  expressible  in  terms  of  the  vectors  a  and  giving 

A  =  52  -  4AC  =  -  (a^S-^a)(/?^S-V)  (6.173) 

The  matrix  is  expressible  as  the  square  of  a  symmetric  real  matrix  in  many,  perhaps  all 
meaningful  cases.  Thus,  we  set 

E"^  = 

The  requirements  for  the  existence  of  the  matrix  f2  are  not  developed  here,  but  assuming 
that  such  an  Q,  is  available  the  vectors 


a  =  Da, 

are  defined.  Then  (6.15)  can  be  written  as 

g  =  D^ 

(6.17) 

rri  rr*  rrt  rr> 

A  =  ga^  g  —  ag  g 

This  is  equivalent  to  a  relationship  between  inner  vector  products 

(6.18) 

II 

< 

-lapIgP 

(6.19) 

which  according  to  the  Schwarz  inequality  cannot  be  positive. 
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FIGURE  6-2.  JOINT  DENSITY  OF  (p,?.r)  ON  A  ^-SECTION,  GENERAL  CASE 


FIGURE  6-3.  AN  ASPECT  OF  THE  JOINT  DENSITY  OF  (p,  q,  r)  ON  A  ^-SECTION,  GENERAL  CASE 
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CHAPTER  7 

THRESHHOLDING 


Let  us  now  consider  threshholding  the  sum-channel  phasors  to  include  only  those  values 
(r,  s)  which  obey  the  condition  r'^  +  >  Rq.  Applying  this  condition  is  cumbersome  unless 

the  density  function  has  a  rather  simple  form.  In  particular,  it  is  straightforward  to  apply  if 
the  means  are  zero.  Consider  first  the  density  function  of  Case  I.  The  threshholding  condition 
is  applied  by  changing  the  limits  on  the  outer  integral  of  (3-5)  from  (0  to  oo)  to  (Rq  to  oo). 
The  result  is 

1  f°° 

- — —  /  R^dRexp 


2^2 


-I- 


2cr2 


(7.1) 


Performing  this  integration  gives 


7^  h?{p^  +  +  7^)  +  1 

TT  (p2  -f  ^2  ^  ^2)2 


exp  +  ?*  +  7^)] 


(7.2) 


where  h  is  written  for  the  quantity  i?o/(<?’xV^)-  This  expression  needs  to  be  normalized, 
because  as  it  stands  it  will  contain  less  than  unit  volume.  The  normalization  can  be  applied 
to  the  modified  density  function  of  the  primary  variables,  or  directly  to  (7.2).  The  latter 
is  more  convenient,  and  accordingly  it  is  found  that  the  integral  of  (7.2)  over  all  p  and  q  is 
exp(— /i^7^).  Dividing  by  this  gives  the  final  result. 


Pp,giP^q) 


7^  +  q‘^  +  j'^)  +  l 

TT  (p2  -f  g2  j2'j2 


exp  {-h^(p^  +  ?^)] 


(7.3) 


Here  also,  the  parameter  7  can  be  hidden  by  appropriate  scaling  of  (4.7)  along  with  the 
definition  h'  =  7/1  to  give 


Pp,g 


(7.4) 


This  radially-symmetric  joint  density  function  is  illustrated  in  Figures  7-1  and  7-2  for  various 
threshhold  values.  Note  that  h'  =  0  corresponds  to  the  case  without  threshholding  for  which 
the  variance  is  infinite.  As  the  threshhold  value  h'  increases,  the  tails  become  shorter  and 
the  variance  decreases.  It  is  rather  complicated  to  obtain  the  marginal  distribution  of  p  or 
q.  However,  it  is  seen  that  the  means  p  and  q  are  zero,  while  the  variance  is  readily  obtained 

,00  .oc 

<^l  =  E[{p-p)^]  =  E\p'^]=  dq  P^Pp,g{p,q)dp  (7.5) 

*/— OO  J —OO 
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Substituting  (7.3)  into  (7.5)  gives  the  variance 

(7.6) 

where  E\  is  the  exponential  integral  defined  by 

7°°  e~^dt 

^l(-2)  =  J  (larg^l  <  tt)  (7.7) 

(see  Abramowitz  and  Stegun  [3]).  This  is  also  the  variance  of  as  can  be  seen  from  the 
symmetry  of  (7.3).  Figure  7-3  shows  the  variance  as  function  of  the  threshhold  value. 


The  same  procedure  is  applied  to  the  integral  of  (5-5)  in  the  case  of  partially-correlated 
isotropic  noise  with  zero  means.  After  the  integration  of  (5.5)  from  0  to  27r  with  respect  to 
dO,  the  unnormalized  result  is 

1  n3  »r.  r  +  7^  n2 

■— ^r— - —  /  WdRexx>\ - t-xt: - - R 

-  pn  Jr 


P?  dR  exp 


'Ro 


27rcr3 
which  is  equal  to 

r hJijP'  q.  ^2  ^  ^2  _  2p7p)  +  1 

(p2  .p  ^2  q.  ^2  _  2/57p)2 

where  h  is  now  given  by 

h 


2cT2(l-p2) 


(7.8) 


21 

IT 


exp  [-h'^ip^  +  9^  +  7^  -  2p7P)] 

Ro 


(Txy/2{1  -  p) 

The  integral  of  (7.9)  over  all  p  and  q  is 

exp[-fi^(l  -  p^)] 

1  —  p2 

The  normalized  density  function  is  (7.9)  divided  by  (7.11),  giving 


(7.9) 


(7.10) 


(7.11) 


PvAp,<i) 


7^(1  -  f^) 


TT 


h?{p^  +  7^  +  7^  -  2/27p)  +  1 

(p2  +  ^2  ^2  _  2p7p)2 


exp[-fi^(p^  -1-  -f  -  2p'yp)] 


2  2 


(7.12) 


or 


7=;, ,(?',!') 


TT 


h'‘2(^pi2  ^  qi2  _  2pp')  -P  1 


exp[-/i'^(p'^  +  9'^  +  -  2pp')]  (7.13) 


(p'2  +  -f  1  —  2pp/)^ 

This  joint  density  function  is  illustrated  in  Figure  7-4.  It  is  seen  by  inspection  that  q  =  0. 
To  obtain  the  mean  value  p  the  substitutions  t  =  p  —  pj  ot  t'  =  tf'j  =  p'  —  p  and  t  =  R  cos  9, 
q  =  Rs\n9  are  made  in  (7.13).  Thus, 


l-p2 


fi'2(t'2  +  qt2^i_p2^1 


K,q  TT  [  (t'2  +  ^'2  +  1  _  ^2)2 


exp[-/i'^(t'^  +  q'^)] 


(7.14) 
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or 


P'  = 
■*  p,g 


l-p2 


TT 


exp(—h'^R^) 


(7.15) 


(R^  +  1-  ^2)2 

This  function  is  shown  in  Figure  7-5.  The  mean  value  t'  is  zero,  and  hence  p  =  p'f,  because 
(7.14)  is  an  even  function  of  t'  and  the  integral  defining  its  mean  converges.  It  may  seem 
somewhat  surprising  that  p  is  independent  of  h. 

Again,  the  marginal  distributions  of  p  and  of  q  are  rather  cumbersome,  and  are  not 
developed  here.  The  variances  of  p  and  of  q,  on  the  other  hand,  are  easily  obtained.  Thus, 


E[t^]  =  1-ii — ^  [  cos^edO  f  R^dR 
Jo  Jo 

which  is  given  by 


h^[R^  +  ^^{l-p^)]  +  l 


[R2+^2(1_^2)]2 


exp[—h^Rj]  (7.16) 


m  = 


7^(1  -P^). 


exp[-/iV(l  -  />^)]jEi[/iV(l  -  P^)] 


2,,2/ 


The  same  result  is  obtained  for  E[q^].  The  variance  of  p  is  then 

W)  =  (^')+pV,  or  E[{p-p^)]  =  E[t^] 


(7.17) 


(7.18) 


The  effect  of  threshholding  can  be  evaluated  by  analysis  only  in  the  simplest  cases.  When 
analytic  solutions  are  not  available  or  are  not  feasible  to  obtain,  Monte  Carlo  methods  can 
be  used  in  the  manner  described  in  Chapter  3.  Rather  severe  threshholding,  with  Ro  =  36, 
was  applied  to  the  example  of  Case  IV  presented  in  Figure  4-10.  The  result  is  shown  in 
Figure  7-6.  The  effect  of  threshholding  is  seen  quite  dramatically  by  comparing  the  two 
figures.  Not  surprisingly,  the  spread  of  points  in  the  threshholding  example  appears  much 
reduced,  and  in  fact,  a  variance  exists  and  can  be  estimated. 

A  similar  study  of  threshholding  with  the  rather  high  value  of  Ro  =  36  was  done  on 
the  example  of  Case  VI  presented  in  Figure  5-8.  The  result,  shown  in  Figure  7-7,  displays  a 
noteworthy  ring  structure.  These  examples  may  be  somewhat  atypical,  as  they  were  selected 
in  order  to  exhibit  such  curious  features.  The  marginal  distribution  of  p  corresponding  to 
Figure  7-7,  obtained  by  330,000  Monte  Carlo  experiments,  is  shown  in  Figure  7-8.  It  is  noted 
that  the  deep  hole,  quite  conspicuous  on  the  bivariate  distribution,  is  much  less  prominent 
on  the  marginal  distribution. 
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FIGURE  7-1.  JOINT  MR  DENSITY  WITH  VARIOUS  THRESHHOLDS 


FIGURE  7-2.  JOINT  MR  DENSITY  WITH  VARIOUS  THRESHHOLDS 
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FIGURE  7-6.  MONTE  CARLO  SCATTER  DIAGRAM  OF  JOINT  MR  DENSITY,  UNCORRELATED  NOISE  WITH  THRESIIHOLD 


FIGURE  7-7.  MONTE  CARLO  SCATTER  DIAGRAM  OF  JOINT  MR  DENSITY  FOR  PARTIALLY  CORRELATED  NOISE  AND  THRESHHOLD 


FIGURE  7-8.  MARGINAL  DISTRIBUTION  OF  REAL(MR)  FOR  PARTIALLY  CORRELATED  NOISE 
AND  TIIRESHIIOLD,  DETERMINED  BY  MONTE  CARLO  METHOD 


NSWCDD/TR-94/97 


CHAPTER  8 
CONCLUSIONS 


The  joint  density  function  of  the  real  and  imaginary  parts  of  the  monopulse  ratio  have 
been  developed  for  various  models  of  the  Gaussian  joint  density  of  the  sum  and  difference 
phasors  of  a  monopulse  radar.  The  resulting  functions  have  been  reduced  to  a  particularly 
simple  form  by  introducing  alternative  parameters  defining  the  joint  density  of  the  primary 
variables  (i.e.,  the  parameters  X,  U,  W,  and  Y  of  Chapter  4).  It  was  shown  that  the 
dependence  on  the  noise  level  in  the  sum  and  difference  phasors  is  usually  representable  in 
terms  of  a  single  parameter,  7,  equal  to  the  ratio  of  the  two  standard  deviations.  A  further 
simplification  in  form  was  achieved  by  scaling  with  factors  dependent  on  7  to  ehminate  this 
parameter  from  the  expressions  defining  the  joint  densities. 

In  a  few  cases  it  was  possible  to  obtain  simple  expressions  for  the  marginal  density  of  the 
real  (or  imaginary)  part  of  the  MR.  Results  were  presented  for  the  cases  of  uncorrelated, 
isotropic  noise  with  zero  means  of  the  primary  phasors  (Case  I),  and  zero  mean  of  the 
difference  phasor  (Case  II),  and  for  partially- correlated  isotropic  noise  with  zero  means  of 
the  primary  phasors  (Case  V).  The  term  “partially-correlated”  is  used  here  to  denote  the 
case  where  the  covariance  matrix  contains  some  zero  elements,  and  does  not  refer  to  the 
degree  of  correlation  or  to  the  value  of  the  correlation  coefficient.  Marginal  densities  of  the 
real  part  of  the  MR  have  been  obtained  in  other  contexts  by  [1,2, 3,4]. 

Generally,  the  marginal  density  functions  of  a  sole  component  of  the  MR  are  more  com¬ 
plicated  and  in  many  cases  less  interesting  than  the  complete  joint  density.  An  expression 
for  the  mean  value  of  the  MR  was  obtained  in  most  of  the  cases.  The  MR  variance,  known  to 
be  infinite  usually  for  an  observation  consisting  of  a  single  pulse,  was  obtained  only  for  the 
cases  resulting  in  a  finite  value,  namely,  where  a  threshhold  on  the  sum-phasor  magnitude 
was  imposed  (in  Chapter  7).  Any  positive  threshhold  value  results  in  a  finite  MR  variance. 

Some  of  the  most  interesting  solutions  are  those  resulting  from  non-vanishing  means  of 
the  primary  variables.  Although  scintillation  and  phase  uncertainty  tend  to  produce  zero- 
mean  density  functions  of  the  primary  variables,  the  consequence  of  non-zero  means  may  be 
applicable  to  some  situations.  One  such  case  is  that  of  observations  on  a  target  located  at  a 
specific  point.  Of  course,  the  analysis  is  much  simpler  if  the  mean  value  of  either  the  sum 
or  the  difference  phasor  is  zero. 

Three  basic  models  of  noise  covariance  were  studied:  uncorrelated  isotropic  noise  (Chapter 
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4),  partially- correlated  isotropic  noise  (Chapter  5),  and  a  general  noise  covariance  matrix 
(Chapter  6).  The  effect  of  the  correlation  coefficient  on  the  joint  MR  density  was  studied. 
Monte  Carlo  simulations  were  used  to  corroborate  the  analytically-obtained  density  functions 
in  a  couple  of  cases.  A  feature  that  has  not  been  previously  noted  in  the  literature  is  a 
prominent  minimum  or  ffiole’  in  the  joint  MR  density  that  occurs  under  some  circumstances. 
Such  a  minimum  has  been  observed  in  cases  of  non-zero  means  of  the  primary  phasors  (as 
in  Figures  4-7  through  4-10  and  5-5  through  5-8).  Some  new  definite  integrals  that  may  be 
useful  in  other  studies  were  evaluated  and  are  given  in  (A.33),  (A.35),  and  (A.36). 

In  carrying  out  this  study  it  became  apparent  to  the  authors  that  a  better  understanding 
of  noise  processes  in  monopulse  radar  systems  is  needed.  The  various  noise  sources  affect 
the  statistical  properties  of  the  density  functions  in  different  ways.  To  our  knowledge,  the 
existing  literature  lacks  a  simple  and  direct  description  of  the  noise  processes.  Development 
of  a  noise  model  which  is  based  on  the  particular  features  of  each  noise  source  is  planned  to 
establish  a  small  number  of  parameters  which  characterize  the  complete  noise  environment 
with  sufficient  accuracy  for  practical  work. 
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APPENDIX  A 


EVALUATION  OF  INTEGRALS 
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The  integration  of  (4.5)  is  carried  out  by  the  substitution  w  =  B?,  which  gives 

fp^  +  q'^  ,  1 


1  f°° 


w 


w  dw 


(Al) 


which  gives  (4.6). 


The  integral  (4.9)  is  evaluated  by  Abramowitz  and  Stegun  ([3],  equation  3.3.24)  who  give 
the  indefinite  integral 


/dx  1  /x\  X 

(rc2  +  a2)  "  lo  J  2a2(a:2  + 


) 


(A2) 


By  setting  +  7^,  the  last  term  vanishes  at  both  limits  yielding  the  result  indicated 

in  (4.9). 

The  integration  of  (4.15)  is  carried  out  by  a  translation  along  the  0  axis  by  expressing  the 
integrand  in  the  form 


exp[C'cos(0  +  a:)]  =  exp[C'(cosQ;cos0  —  sinasin^)] 

where  C  and  a  are  to  be  determined.  Evidently  they  must  satisfy  the  relations 

Ccosa  =  Ru  I 
Csina  =  —Rvla\ 

which  determine  the  value  of  C  to  be 

^  Ry/u^  + 

^  “  2 

A  transformation  is  now  made  to  a  new  variable  of  integration, 

^  =  6  a 


(X.3) 


(AA) 


(A.i) 


(A6) 


The  angle  a  can  similarly  be  determined,  but  its  value  is  irrelevant  as  the  integration  is  over 
a  complete  cycle  of  either  6  or  </>.  The  result  is 

r2w 

(A7) 


I  exp[C  cos(l>]d(f> 
Jo 


which  has  the  value 


27r/o(C')  =  2wIo 

as  given  by  [3]  (equation  9.6.19). 


Ry/ 1/2  -|- 1;2 


(.4,8) 


The  authors  are  indebted  to  John  E.  Gray  for  giving  direction  in  the  evaluation  of  this 
integral.  The  integral  of  (4.17)  is  carried  out  by  expressing  the  integral  in  the  form 


r°°  9 

/  R^e-^^  h{BR)dR 

Jo 


(A.9) 
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where  A  and  B  are  constants  with  respect  to  the  integration,  and  are  readily  evaluated  in 
terms  of  the  parameters  in  (4.17).  Equation  (9.6.10)  of  [3]  gives  the  following  series  expansion 
for  the  Modified  Bessel  Function, 


k=0 


(.4.10) 


Substitution  into  (A. 9)  gives 


“2-2^52^  roo 

Jo 


(A.ll) 


The  integral  in  the  above  expression  can  be  evaluated  by  making  the  substitution  2:  =  AR?, 
which  gives 

The  integral  part  of  (A. 12)  is  just  (k  +  1)!,  as  can  be  found  in  any  standard  integral  table. 
Substituting  back  into  (A.ll)  gives 


^2-2*=52^(jb+l)!  1 

2-j  9W\]2Ak+2  ~0A2l^ 


k-=0 


2[it!]2A^+2  2A2  ^  k\ 


4A, 


1 


2A2 

1 


1  f  B^\^~\  ^  1_  f 


/ B^\  i  [  ^  \ 

Wj  (TH)!  {UJ  +  ^  W  K4AJ 


(A13) 


52 

exp  I  7T 


(^)  A  1 


iAj  LV4>17 


Considering  (4.18)  to  show  that  the  integral  over  all  p  and  q  is  unity,  first  make  the 
substitution 

p  =  R  cos  9 

R  -  , 

q  =  Rsm  9 


Thus, 


/OO  ^00  fOO  r2lt  roo 

dq  Pp,q{p^q)dp=  RdR  Pp,g{p,  q)d9  =  2Tr  RdRPp^g{p,q)  (A.15) 
-OO  J —OO  Jo  Jo  Jo 

To  facilitate  the  integration  with  respect  to  R,  p  is  used  as  the  variable  of  integration,  where 

7/72 


p  = 


The  integral  then  becomes 


e-^ 

^  (1  +  p)  e^dp 

A-4 


(A.16) 


(A.17) 
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which  is  unity. 


To  evaluate  the  integral  (4.22)  the  exponential  is  expressed  as  a  power  series  and  the 
integration  is  performed  term  by  term.  For  convenience  a  change  of  variable  is  made  so  that 
the  quadratic  term  in  the  radical  will  have  a  coefficient  of  unity.  Thus, 


z  =  ny/p^  +  7^ 

Evaluation  of  each  term  in  the  resulting  series  requires  the  integral 

z'^dz 


i: 


(A18) 


(A19) 


z  —  Z^ 


for  n  =  1,2, ....  The  integral  exists  even  though  the  denominator  in  the  integrand  vanishes 
at  both  limits.  The  result  is 

7r(2n-l)U» 

2«-i(n-l)!n!  ^  ‘  ^ 

When  the  exponential  factor  is  included,  the  result  can  be  expressed  in  series  form.  Thus, 
r2A 


I 


oo  r2A 


0  V2Az  —  z"^  Jo  V2Az  — 


^  fAY^  (2n 


+  2k-l)\{ABy 


(A.21) 


2”+^/:!(n  +  k  —  l)!(n  +  fc)! 

With  this  result  the  integral  (4.22)  can  be  readily  evaluated. 

To  show  that  the  density  function  Pp{p)  of  (4.23)  encloses  unit  area,  the  integral 

dp 


/: 


1-00  (p2  -)-  ^2^n+2 

is  needed.  With  the  substitution 

p  =  7  tan  9,  dp  =  'f  sec^9  dO 

the  integral  can  be  written 

7  sec^O  dO 


(A.22) 


(A.23) 


/•7r/2 

'I 


72n+l(tan^0  +  1)"+ 


/•7r/2 

-  =  27-2”  /  cos2”-i 

2  Jo 


Ode 


(A.24) 


The  definite  integral  occurring  in  (A.24)  is  found  in  any  standard  table  of  definite  integrals, 
with  the  final  result  being 


27 


— 2n 


2-4-6---(2n-2) 


[l-3-5---(2n-l)J 

A-5 


22n-l[(„_l)!]27-2n 

(2n-  1)! 


(A.25) 
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The  expression  (4.23)  can  be  written 


dp 


J=0  \  /  ./-CX)  (p2  -f  ^2p-|-l+2 

When  (A.25)  is  substituted  into  (A.26)  the  integral  can  be  written 


(A.26) 


/OO  °°  Tjj 

Pp{p)dp  =  TT  =  ^ 

jk=o 


(A.27) 


An  alternative  formulation  in  terms  of  standard  functions  is  obtained  with  the  help  of  formula 
No.  3.364.1  of  [4],  which  states  that 


I 


^  e-P^dx 
0  y/x{2  -  x) 


=  Tre  ^Io{p)  for  p  >  0 


(A.28) 


It  can  be  demonstrated  that  the  restriction  on  p  can  be  relaxed  to  include  all  real  values. 
With  the  substitution  p  —  2  —  x  the  integral  (A.28)  becomes 


-2p 


I 


or 


I 


2  eP'^dp 


=  Tre  ^/o(p)  for  p  >  0 


=  Tre^/o(p)  for  p  >  0 


(A.29) 


(A.30) 


/o 

Since  /o  is  an  even  function  this  result  essentially  shows  that  (A.28)  is  valid  for  negative  as 
well  as  positive  p,  while  the  value  p  =  0  can  obviously  be  included  for  continuity  reasons. 
(Skeptics  can  note  that  p  =  0  gives  a  standard  integral  found  in  elementary  tables,  which  is 
consistent  with  (A.28)  and  (A. 30).) 


The  next  step  is  to  further  generalize  (A.28)  by  making  the  substitutions 


_  2{ii  -  A)  _  -a{B  -  A) 

B-A  2 


where  A  <  whereupon  the  integral  becomes 


—aA 


e°‘>^du  { 

/  — ,  . . =  TT  exp 

Ja  \ 


o(B  -  A)\  /  a(B-A) 


or 


e°‘^da  fa(A  +  B)\  f  a{B  —  A) 

L  J(,-A)(B-,)  =  "  ( - r- 


(A31) 


{A.Z2) 


{A.33) 


Integrals  containing  powers  of  p  in  the  integrand  can  be  obtained  from  (A.33)  by  differenti¬ 
ation  with  respect  to  a.  Derivatives  of  the  Modified  Bessel  Functions  are 


l'o{z)  =  hiz),  I[{z)  =  Io{z)- 


h{z) 


(A.34) 
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which  give  the  result 


f 


A  ^/if^  -  A){B  -  n) 


fa{A  +  B)\ 

=  TT  exp  - - -  1  X 


A  +  B\  /  a{B-Ay 


(^.35) 


and 


f 


=  Trexp 


a{A  +  B) 


)l( 


U  -  A){B  -  fi)  V  2 


A^  +  B^\  ,  f-a{B-A)' 
h 


(A.36) 


These  formulas  are  used  to  obtain  the  expression  (4.25)  by  setting  A  =  0,  B  =  lijp^  +7^)) 
and  a  =  1. 

To  show  that  the  joint  density  function  (4.32)  contains  unit  volume,  first  make  the  trans¬ 
formation  (A. 14),  whereupon  the  integral  over  all  p  and  q  becomes 

/oo  yoo  yoo  f2ir  yoo 

dp  Pp,qiP^q)dq=  RdR  Pp^g{p,q)d9  =  2ir  Pp^g{p,q)RdR  (A37) 
-00  J —<x>  Jo  Jo  Jo 


Next,  transform  the  variable  of  integration  in  (4.32)  by  letting 


a  =  p  —  X  = 


^2+^2 


and  the  integral  becomes 


(l-hX-f-a)e“da 

A  J-X 

which  is  readily  integrated  and  is  equal  to  unity,  as  required. 


(A.38) 


(A.39) 


The  normalization  of  the  expression  (7.2)  is  carried  out  by  the  transformation  (A. 14), 
which  gives 


poo  poo  pc 

/  /  Pp,g{p^q)^pdq  =  ^i^ 

J—oo  J — oo  •'0 


With  the  substitution 


this  becomes 


fe2(^2  ^  ^2)  1 

722+^2 


exp[-;i2(;^2  ^  ^2)1^  (^  40) 


t  = 


(A41) 

(A.42) 
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These  integrals  are  expressible  by  the  exponential  integral  functions  described  in  [3],  which 
in  terms  of  (A.42)  becomes 

^^h^E,{'r^h‘^)  +  E2ij^h^)  (A43) 

There  is  the  recursion  relationship 

E2{z)  =  e-^-zEi{z)  (A.44) 


which  allows  the  E2  function  of  (A.43)  to  be  expressed  in  terms  of  the  Ei  function,  which 
cancels  with  the  first  term  and  gives  the  result  exp(— as  indicated  in  Chapter  7. 

The  variance  is  given  by  (7.6),  which  is  transformed  by  the  substitution  (A. 14)  giving 


(722  +  ^2)2 


exp{-h‘^R^)dR 


With  the  substitution  (A.41)  this  becomes 


(A.45) 


+  (1  -  (A.46) 

2  Ji 

The  second  and  third  terms  of  (A.46)  are  again  evaluated  in  terms  of  the  exponential  integral 
functions  to  obtain 

o-p  =  +  (1  -  h^^‘‘)Ei(hW)  -  E2(hW)] 

And  again  the  E2  function  is  expressed  in  terms  of  the  E\  function  by  (A.44)  to  give  the 
result  indicated  in  (7.7).  The  same  steps  are  used  in  the  evaluation  of  (7.20). 


Integration  of  (6.7)  is  facilitated  by  the  substitutions  of  (6.9)  giving 


^P,g,r(p,g,0  = 


— G  roo 


(A.48) 


with  f  as  defined  in  (6.9).  Now,  change  of  integration  variable  from  stow  via  the  substitution 


5  +  r  =  ly 


(A.49) 


gives 

-G  foo  „ 

-  2ru,  +  (f^  +  dy,  (.4.50) 

The  middle  term  vanishes  because  it  is  an  odd  function.  The  result,  (6.8),  is  obtained  from 
any  standard  table  of  definite  integrals.  Then  representing  G  in  powers  of  r, 
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defines  the  new  parameters, 


A  =  52  -  4AC 


H 


A 

4(7 


M  =  F-HL^-^  =  F  +  A-\AF^  -  BDF  +  CD^) 


(A.52) 


The  marginal  density  Pp^q  then  can  be  obtained  by  integrating  the  expression  (6.8)  with 
respect  to  r  to  obtain 


„  ,  ,  (7r(7)-3/2e-^  f°°  r/R2  \ 


EB 


8\/i2l 

(7r(7)-3/2e-^  r/^ 

7-00  LV2C' 


r  +  1:^  +  1 


C 


F^ 


8v1q 

where  a  change  of  variable 


+  2C^  v?  +  N 


w  =  r  A  L 


2(7 


e-^^^dw 


exp[— i7(r  +  L)^]dr 

(A.53) 

(A.54) 


has  been  made.  (The  term  that  is  linear  in  w  has  been  omitted  since  it  contributes  nothing 
to  the  integral.)  The  result  is 


Pp,q{P^<l)  = 


(iVA  -B^-  4C2)e-^ 

4C7^A^/^^\^A 


(A.55) 


where 


N  =  A-^[2B^CD^  +  2B^CF^  -  SBC^DF  +  SC^D^  -  8ABCDF 
A  SA^CF^  +  -  SAB^C  A  16A^C^] 


(AM) 


with  A  as  defined  in  (6.12)  and  (6.14).  Finally,  (A. 55)  reduces  to  (6.11). 
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